Implementation of on-site velocity boundary conditions for D3Q19 lattice Boltzmann

On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.Comment: 13 pages, 4 figures, revised versio

Interface Roughening in a Hydrodynamic Lattice-Gas Model with Surfactant

Using a hydrodynamic lattice-gas model, we study interface growth in a binary fluid with various concentrations of surfactant. We find that the interface is smoothed by small concentrations of surfactant, while microemulsion droplets form for large surfactant concentrations. To assist in determining the stability limits of the interface, we calculate the change in the roughness and growth exponents $\alpha$ and $\beta$ as a function of surfactant concentration along the interface.Comment: 4 pages with 4 embedded ps figures. Requires psfig.tex. Will appear in PRL 14 Oct 199

Towards Better Integrators for Dissipative Particle Dynamics Simulations

Coarse-grained models that preserve hydrodynamics provide a natural approach to study collective properties of soft-matter systems. Here, we demonstrate that commonly used integration schemes in dissipative particle dynamics give rise to pronounced artifacts in physical quantities such as the compressibility and the diffusion coefficient. We assess the quality of these integration schemes, including variants based on a recently suggested self-consistent approach, and examine their relative performance. Implications of integrator-induced effects are discussed.Comment: 4 pages, 3 figures, 2 tables, accepted for publication in Phys. Rev. E (Rapid Communication), tentative publication issue: 01 Dec 200

Multi-component lattice-Boltzmann model with interparticle interaction

A previously proposed [X. Shan and H. Chen, Phys. Rev. E {\bf 47}, 1815, (1993)] lattice Boltzmann model for simulating fluids with multiple components and interparticle forces is described in detail. Macroscopic equations governing the motion of each component are derived by using Chapman-Enskog method. The mutual diffusivity in a binary mixture is calculated analytically and confirmed by numerical simulation. The diffusivity is generally a function of the concentrations of the two components but independent of the fluid velocity so that the diffusion is Galilean invariant. The analytically calculated shear kinematic viscosity of this model is also confirmed numerically.Comment: 18 pages, compressed and uuencoded postscript fil

Pipe network model for scaling of dynamic interfaces in porous media

We present a numerical study on the dynamics of imbibition fronts in porous media using a pipe network model. This model quantitatively reproduces the anomalous scaling behavior found in imbibition experiments [Phys. Rev. E {\bf 52}, 5166 (1995)]. Using simple scaling arguments, we derive a new identity among the scaling exponents in agreement with the experimental results.Comment: 13 pages, 3 figures, REVTeX, to appear in Phys. Rev. Let

Fluctuations of elastic interfaces in fluids: Theory and simulation

We study the dynamics of elastic interfaces-membranes-immersed in thermally excited fluids. The work contains three components: the development of a numerical method, a purely theoretical approach, and numerical simulation. In developing a numerical method, we first discuss the dynamical coupling between the interface and the surrounding fluids. An argument is then presented that generalizes the single-relaxation time lattice-Boltzmann method for the simulation of hydrodynamic interfaces to include the elastic properties of the boundary. The implementation of the new method is outlined and it is tested by simulating the static behavior of spherical bubbles and the dynamics of bending waves. By means of the fluctuation-dissipation theorem we recover analytically the equilibrium frequency power spectrum of thermally fluctuating membranes and the correlation function of the excitations. Also, the non-equilibrium scaling properties of the membrane roughening are deduced, leading us to formulate a scaling law describing the interface growth, W^2(L,T)=L^3 g[t/L^(5/2)], where W, L and T are the width of the interface, the linear size of the system and the temperature respectively, and g is a scaling function. Finally, the phenomenology of thermally fluctuating membranes is simulated and the frequency power spectrum is recovered, confirming the decay of the correlation function of the fluctuations. As a further numerical study of fluctuating elastic interfaces, the non-equilibrium regime is reproduced by initializing the system as an interface immersed in thermally pre-excited fluids.Comment: 15 pages, 11 figure

Casimir energy of massive MIT fermions in a Bohm-Aharonov background

We study the effect of a background flux string on the vacuum energy of massive Dirac fermions in 2+1 dimensions confined to a finite spatial region through MIT boundary conditions. We treat two admissible self-adjoint extensions of the Hamiltonian and compare the results. In particular, for one of these extensions, the Casimir energy turns out to be discontinuous at integer values of the flux.Comment: 16 pages, 3 figure

Massive 3+1 Aharonov-Bohm fermions in an MIT cylinder

We study the effect of a background flux string on the vacuum energy of massive Dirac fermions in 3+1 dimensions confined to a finite spatial region through MIT boundary conditions. We treat two admissible self-adjoint extensions of the Hamiltonian. The external sector is also studied and unambiguous results for the Casimir energy of massive fermions in the whole space are obtained.Comment: 12 pages, 5 figures, LaTe

Spinodal decomposition of off-critical quenches with a viscous phase using dissipative particle dynamics in two and three spatial dimensions

We investigate the domain growth and phase separation of hydrodynamically-correct binary immiscible fluids of differing viscosity as a function of minority phase concentration in both two and three spatial dimensions using dissipative particle dynamics. We also examine the behavior of equal-viscosity fluids and compare our results to similar lattice-gas simulations in two dimensions.Comment: 34 pages (11 figures); accepted for publication in Phys. Rev.

Density waves and 1/f1/f density fluctuations in granular flow

We simulate the granular flow in a narrow pipe with a lattice-gas automaton model. We find that the density in the system is characterized by two features. One is that spontaneous density waves propagate through the system with well-defined shapes and velocities. The other is that density waves are so distributed to make the power spectra of density fluctuations as $1/f^{\alpha}$ noise. Three important parameters make these features observable and they are energy dissipation, average density and the rougness of the pipe walls.Comment: Latex (with ps files appended